Proximity-driven source of highly spin-polarized ac current on the basis of 
superconductor/weak ferromagnet/superconductor voltage-biased Josephson junction. 
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We theoretically investigate an opportunity to implement a source of highly spin-polarized ac cur- 
rent on the basis of superconductor/weak ferromagnet/superconductor (SFS) voltage-biased junc- 
tion in the regime of essential proximity effect and calculate the current flowing through the probe 
electrode tunnel coupled to the ferromagnetic interlayer region. It is shown that while the polariza- 
tion of the dc current component is generally small in case of weak exchange field of the ferromagnet, 
there is an ac component of the current in the system. This ac current is highly spin-polarized and 
entirely originated from the non-equilibrium proximity effect in the interlayer. The frequency of the 
current is controlled by the voltage applied to SFS junction. We discuss a possibility to obtain a 
source of coherent ac currents with a certain phase shift between them by tunnel coupling two probe 
electrodes at different locations of the interlayer region. 

PACS numbers; 74.50.-|-r, JiAb.+c 
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Various proximity and transport phenomena in hybrid 
structures containing superconducting and ferromagnetic 
elements are currently in the spotlight. The equilib- 
rium transport and proximity effect in such structures 
have been theoretically and experimentally investigated 
recently in details as for the case of weak ferromagnetic 
alloys (see Refji and references therein) so as for half- 
metals like CrO^^. 

Spin-dependent properties of nonequilibrium systems 
are also actively investigated. In particular, the spin im- 
balance induced in a ferromagnet (F) in nonequilibrium 
conditions was studied in ferromagnet-superconductor- 
ferromagnet (FSF) junctions^i^iii^i^. It was found that 
antiferromagnetic alignment of the exchange fields of 
the ferromagnets strongly suppresses superconductivity, 
which leads to large magnetoresistive effect. The influ- 
ence of the interplay between Andreev reflection at the 
interface and spin accumulation close to the interface was 
investigatediS. The effect of spin injection in Josephson 
junctions was considered^^. The papers^' have theoret- 
ically studied the spin relaxation due to spin-flip scatter- 
ing. 

The authors o f^'^d^ propose the possibility of manipu- 
lating magnetization of a mesoscopic normal (N) region 
through the Zeeman splitting of superconducting den- 
sity of states and applied voltage in voltage biased SNS 
or FSNSF tunnel junctions. Also, the spin-polarized 
transport through superconductor-normal metal hybrid 
structures, where the density of states in a supercon- 
ductor is Zeeman-splitted, was investigatedi^iiiii^. It 
was shown that, in alternative to half- metallic ferro- 
magnets, such junction can be used to generate highly 
spin-polarized currents, which are tunable in magni- 
tude and sign by the bias voltage and exchange field. 
j ^i3,i4 ^]^g jjormal metal region between the supercon- 
ducting leads has been considered to be long enough 
in order to suppress the proximity effect and, conse- 
quently, ac Josephson effect in it. However, if the junc- 
tion length is of the order of superconducting coherence 



length ^ = ^hD/A, where D is the diffusion constant 
and A is the superconducting order parameter in the 
leads, the interplay between the proximity effect and 
the Zeeman exchange field in the interlayer can lead to 
a number of novel qualitative effects. In particular, it 
was shown that by tunnel coupling of the normal re- 
gions of two superconductor-normal metal-ferromagnet 
trilayers the absolute spin- valve effect can be realized for 
a certain interval of voltages applied between the nor- 
mal regions^-. For voltage-biased superconductor-weak 
ferromagnet-superconductor Josephson junction the in- 
terplay between the proximity effect and the Zeeman ex- 
change field results in additional peak-like features in 
the I-V characteristics^^. The ac Josephson current in 
superconductor-ferromagnet-superconductor voltage bi- 
ased junction was also studied^^. 

In the present paper we study voltage-biased 
superconductor-weak ferromagnet-superconductor (SFS) 
Josephson junction focusing on the short enough inter- 
layer (which is considered to be of the order of super- 
conducting coherence length). In this regime the prox- 
imity effect between ferromagnet and superconductors is 
essential. One of the well-known manifestations of the 
proximity effect, which already takes place in equilib- 
rium, is the so-called minigap in the local density of states 
(LDOS) of the interlayeriSiSfi,. Another important man- 
ifestation of the proximity effect is an ac current, which 
appears in case of voltage-biased junction. We show that 
the interplay of two above mentioned phenomena in non- 
equilibrium SFS Josephson junction gives an opportu- 
nity to implement a source of highly spin-polarized ac 
current by tunnel coupling the interlayer region to an 
additional electrode. The frequency of this ac current is 
controlled by the voltage V, applied to the SFS junction. 
In addition the proximity effect in the interlayer causes 
substantial non-linearities in the I — Vp characteristics of 
the ac current flowing through the additional electrode. 
Here Vp is the potential of the probe electrode. The dc 
current flowing through the additional electrode is also 
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considered, although its spin polarization is obtained to 
be rather weak except for some narrow ranges of Vp. We 
also discuss the phase difference between the ac current 
flowing through the probe electrode and the ac Josephson 
current flowing across the SFS junction and show that it 
depends on the position in the interlayer. Therefore the 
system under consideration can be used as a source of 
phase-shifted coherent ac currents by tunnel coupling of 
two probe electrodes at different locations of the inter- 
layer region. 

Further the model and the method we use are de- 
scribed. We study a voltage-biased SFS junction, where 
F is a diffusive weak ferromagnet of length d coupled 
to two identical superconducting reservoirs. The su- 
perconductors are supposed to be diffusive and have s- 
wave pairing. We assume the SF interfaces to be not 
fully transparent and suppose that the resistance of the 
SF boundary Rg dominates the resistance of the fer- 
romagnetic interlayer Rp. We assume the parameter 
(Rp / Rg){o'F / <ys), where ap and Cs stand for conduc- 
tivities of ferromagnetic and superconducting materials 
respectively, to be also small, what allows us to neglect 
the suppression of the superconducting order parameter 
in the S leads near the interface. In addition, a normal 
voltage-biased "probe" terminal is tunnel coupled to the 
interlayer through a junction of a resistance Rp ^ Rg. 

We use the quasiclassical theory of superconductiv- 
ity for diffusive systems in terms of time-dependent Us- 
adel equations^i. The fundamental quantity for diffusive 
transport is the momentum average of the quasiclassi- 
cal Green's function g{x,e,t) = {g{pf,x,£,t))pj. It is 
a 8 X 8 matrix form in the product space of Keldysh, 
particle- hole and spin variables. In general the quasi- 
classical Green's functions depend on space -R, time t 
variables and the excitation energy e. The considered 
problem is effectively one-dimensional and R = x, where 
cc - is the coordinate measured along the normal to the 
junction. 

The quasiclassical Green's function g{x, e, t) satisfies 
the non-stationary Usadel equation, which in the ferro- 
magnetic region takes the form 

[eT3-t{x,e,t),g]^ + ^d,ig®do,g) =0 , (1) 
supplemented by the normalization condition 
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The product ® of two functions of energy and time is 
defined by the noncommutativc convolution A (Ei B — 
e'{dtdt-d^S^)/-2A{£,t)B{e,t). In the problem we con- 
sider the self-energy takes the form S(a;,e,f) = ha^, 
where h is an exchange field in the ferromagnet. and 
di are Pauli matrices in particle-hole and spin spaces, 
respectively. 

In order to solve the Usadel equation it is convenient to 
express quasiclassical Green's function g in terms of Ric- 
cati coherence functions ^^-'^ and which measure 



the relative amplitudes for normal-state quasiparticle and 
quasihole excitations and distribution functions x^ and 
x^ . All these functions are 2x2 matrices in spin space 
and depend on [x,e,t). The corresponding expression for 
g takes the form^^ 
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Riccati coherence and distribution functions obey 
Riccati-type transport equation a^^'^'^ . For the considered 
problem in the ferromagnetic region the equations read 
as follows 

2£®7^- [/i(T3,7^]- 
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Particle-hole conjugation, denoted by ~ , is defined by the 
operation a{x^e^t) — a{x, —e,t)* . In addition to the 
conjugation symmetry, the coherence and distribution 
functions obey the following symmetries 7'^(x,e,t) = 
7^(x, e, ty, x^{x, e, t) = x^{x, e, . 

Eqs. (H]),® should be solved together with the bound- 
ary conditions at SF interfaces. As it was mentioned 
above we consider the case when the dimensionless con- 
ductance of the boundary G — Rp/Rg < 1, so the 
interface transparency Tsp ^ G{l/d) 1. Due to 
the smallness of the interface transparency T we can 
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use Kupriyanov-Lukichev boundary conditions at SF 
boundaries^. In terms of Riccati coherence and distri- 
bution functions they take the form 

Oxll,r [iS;Lr + ll,r ® 9S;Lr- 



<9S;Lr+lLr 



(8) 
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Here Riccati coherence and distribution functions de- 
noted by the lower case symbols I, r are taken at the left 
and right ends of the ferromagnet. The quantities de- 
noted by the lower case symbols {S; I, r) are correspond- 
ing Green's functions at the superconducting side of the 
left and right SF interfaces. As it was already mentioned 
above, under the condition [Rp / Rg){(7 p / a s) ^ 1 we can 
neglect the suppression of superconducting order param- 
eter in the leads and, moreover, take the Green's func- 
tions at the superconducting side of the boundaries to be 
equal to their bulk values, which can be easily deduced 
from the expressions ([3]) and ([4]) using the following bulk 
values of Riccati coherence and distribution functions 
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use of Eq. ((T2)) with the substitution (T3 for ctq. = 1/2 
is the electron spin. In Eq. Vf^y is the Fermi velocity 
component normal to the junction between the interlayer 
and the probe electrode, Nf is the density of the states on 
the Fermi level in the probe electrode, (...)p^ = / dVl/An 

FS 

denotes average over the Fermi surface. 

In the first order on the transparency Tpp of the 
junction between the interlayer and the probe electrode 
(which is considered to be small) the difference between 
the Green's functions for an incoming pj and outgoing 
Pj, quasiparticle trajectories gp{pf,x, e, t) — gp{p pX, e, t) 

entering Eq. \12\ can be expressed in terms of the Green's 
functions corresponding to the uncoupled interlayer and 
probe electrode regions as follows^ - 
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Here the Green's function gp{pf, x, e, t) in the dirty fer- 
romagnetic interlayer is only slightly dependent on the 
momentum direction on the Fermi surface and approxi- 
mately equal to its momentum average value gp{x, e,t). 
[A, B]^ means the commutator A^B — B^A. Substitut- 
ing Eq. p3|) into the expression for the current and 
taking into accout the explicit form of the Green's func- 
tion for the uncoupled normal probe electrode g^^^'^ = 



^piTrfsCTo and gp 
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-2z7rtanh[(£ — eVp)/2T]T3ao, we 



get the current flowing through this electrode: 
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A is the superconducting order parameter absolute value 
in the bulk, which is assumed to be the same in the both 
superconductors. V; is the electric potential in the bulk 
of left (right) superconductor, soV = Vr — Vi is the volt- 
age bias applied to the junction. T is the temperature. 

The electric and spin currents flowing through the 
probe electrode should be found via Keldysh part of the 
quasiclassical Green's function. It is convenient to cal- 
culate the currents in the probe electrode. Given the 
quasiclassical Green's function gp{pf,e,t) at the inter- 
face of the probe electrode the corresponding expression 
for the electric current reads as follows 
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where e is the electron charge and h = \ throughout the 
paper. The spin current j'^^/s'^ can be calculated making 



tanh ■ 
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) gp{x, £, t) — gp{x, £, t) (g) tanh ■ 



2T 



(14) 

cjp''^'^ {x,e,t) is the upper left part of the interlayer 
Green's function cjp'^'^ (x, e, t) in the particle-hole space. 



The spin current j^"^ j s'^ can be calculated from Eq. ([M]) 
with the substitution (73 for iTq . In Eq. (|14p the resistance 
of the interface between the interlayer and the probe elec- 
trode Rp^ = l^,^^^^[d^lAi:)vf,yNj,yTpp. 

The Green's functions g^"^"^ (x,£,t) are calculated 
making use of the Riccati-parameterization technique de- 
scribed above. The time dependence of the supercon- 
ducting order parameter in the leads ~ Ae^**^^*, which 
cannot be removed by the gauge transformation in the 
case of voltage-biased SFS junction, give rise to time de- 
pendence of the Green's function in the interlayer: 
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If the interlayer is long in comparison with the coherence 
length, all the Green's function harmonics corresponding 
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to m 7^ are negligible. This leads, in particular, to the 
suppression of the ac Josephson effect in SFS junction. 
However, if the length of the interlayer is of the order of 
the coherence length, the non-zero harmonics are impor- 
tant and give rise not only to ac Josephson effect, but 
also to ac electric and spin currents flowing through the 
probe electrode. It is worth to note here that we con- 
sider exchange fields, which are very weak as compared 
to the Fermi energy sp and calculate all the currents to 
zero order of the parameter h/ep. In this approxima- 
tion the spin current across SFS junction is absent (it 
only appears in the first order of this parameter). At 
the same time non-zero spin current flows through the 
probe electrode. The qualitative reason for this effect 
is the manifestation of the proximity induced minigap, 
which is splitted by the exchange field, in the essentially 
nonequilibrium electron distribution function in the in- 
terlayer. This is discussed in detail below. 

At first, let us focus on the behavior of the dc compo- 
nent of the electric and spin currents flowing through the 
probe electrode. The differential conductance dl/dV of 
the electric current (normalized to its asymptotic value 
l/Rp corresponding to high enough \Vp\) in dependence 
on Vp is plotted in Figs. [Ija) and ^a). The offset is 
for clarity. Figs.[Hb) and [2Kb) represent the dependence 
of the dc spin current on Vp. Fig. [2] demonstrates how 
these quantities are affected by the exchange fleld, while 
Figs. [Ua)-(b) show the influence of voltage V applied to 
SFS junction. 

The most interesting feature of the electric current is 
series of dips in the differential conductance, which is 
seen in Figs.[lja) and[2Ka). These dips are direct conse- 
quence of the proximity effect and reflect the spin-split 
minigaps in the LDOS extended from the left interface 
(corresponding to the voltages Vp = ±/i) and from the 
right one (located at the positions Vp — V zt h). The 
sign ± is related to the different spin subbands. It is 
seen in Fig. [TJa) that the relative positions of the dip 
pairs can be adjusted by manipulating the voltage V ap- 
plied to SFS, which is easily controlled experimentally. 
Fig. ^a,) demonstrates that these features are the most 
pronounced for weak exchange fields h < A. For h > A 
the minigaps are pushed out from the subgap regions 
of the LDOS and convert into obscure features. Conse- 
quently, the dips in the differential conductance become 
evanescent for high enough exchange fields. 

For weak exchange fields h < A the dc spin current, 
which is represented in Figs. [IJb) and (Hb), is also a 
highly non-linear function of eVp if it is roughly between 
—A and A -|- eV. This is again the consequence of the 
interplay between the minigaps extended from the both 
interfaces of the ferromagnet. It is worth to note here 
that the whole picture is symmetrical with respect to 
Vp — V/2 ^ —{Vp — V/2). The reason is that the curves, 
represented in Figs. [1] and [2l are obtained at a: = d/2, 
that is exactly in the middle of the interlayer. If one 
would calculate the currents flowing through the probe 
electrode closer to one of the interfaces, there would be 
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FIG. 1: (a) The differential conductance of the dc electric 
current flowing through the probe electrode as a function of 
eVp. Different curves are related to three difi^erent values of 
the voltage V, applied to SFS junction. The conductance 
is normalized to its asymptotic value 1/Rp corresponding to 
high enough \Vp\. The offset is for clarity, (b) The depen- 
dence of the dc spin current on eVp for the same three values 
of V. There is no offset. The spin current is measured in 
units of Tc/Rp, where Tc is the critical temperature of the 
superconducting leads, (c) The degree of spin polarization 
P''*''-' of the dc current for eV = 2.5Tc. For other values of the 
voltage V the behavior of p('*'^' is qualitatively the same. The 
other parameters of the junction are the following: h = O.QTc, 
d = 0.94^, G = 0.04. All the results presented in the paper 
are calculated for the temperature T = O.OlTc. 



an asymmetry of the corresponding curves with respect 
to Vp — V/2 — » —{Vp — V/2). This asymmetry origi- 
nates from the fact that the proximity features extended 
from the nearest boundary dominate the proximity fea- 
tures, which extended from the other one. It is seen in 
Fig. [IKb) that the distance between the non-linear fea- 
tures in the spin current is again can be controlled by 
the voltage V. Another non-trivial characteristic feature 
is that the dc spin current tends to a constant value at 
large enough | Vp \ instead of to be a linear function of the 
voltage bias. This constant value is a non-monotonous 
function of the exchange field and declines upon increas- 
ing h (just as the amplitude of the non-linear features 
does). This behavior is a consequence of the manifesta- 
tion of spin-split minigap in the nonequilibrium quasi- 
particle distribution function in the interlayer and is dis- 
cussed qualitatively below. However, it is obvious that 
full dc spin current does not vanish upon increasing h if 
one takes into account the contributions of the first and 
the following orders of h/ep, which are disregarded here 
but become essential for larger exchange fields. 

Fig. [ijc) demonstrates the degree of spin polarization 
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FIG. 2: (a) The differential conductance of the dc electric 
current flowing through the probe electrode as a function of 
eVp. Different curves are related to four different values of 
exchange field h. The offset is for clarity, (b) The dependence 
of the dc spin current on eVp for the same values of h. There 
is no offset. The other parameters of the junction are the 
following: eV = 2.5Tc, d = 0.94^ and G = 0.04. 



of the current as a function of Vp. For the dc component 
of the current this quantity is defined as 



p(dc) 



I JT - ii 



(16) 



where is dc component of the current of spin-up 

(spin-down) electrons, which can be calculated according 
to Eq. dH]) with the substitution {ao+a^j/i (((7o-(73)/2) 
for (To . It is seen that the polarization of the dc current is 
weak everywhere except for the particular narrow ranges 
of Vp, where the electric current is small due to smallness 
of the voltage bias Vp — V/2 applied between the interlayer 
and the probe electrode. 

Now we turn to the discussion of ac current. Fig.[3]is a 
representative example of the dependencies of ac electric 
and spin currents on Vp. Panel (a) demonstrates the cur- 
rents in the middle of the interlayer (at x — d/2), while 
the currents at the left interface {x = 0) are plotted in 
panel (b). We only consider the first harmonics of ac 
current, corresponding to m = ±1 in Eq. psp . The fol- 
lowing harmonics are negligible because we assume the 
transparency of SF interfaces to be low enough, what 
leads to the suppression of a Green's function compo- 
nent gm+i by a factor of G with respect to gm- There- 
fore, the electric and spin ac currents can be expressed 
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Then the 



amplitudes of these currents, which are represented in 
Fig. [3 take the form | = ^^^"^^1 + 
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FIG. 3: The dependencies of the ac electric (solid line) and 
ac spin (dashed line) currents on Vp. h — 0.9Tc, V = 2.5Tc, 
other parameters of SFS junction are the same as in Figs. [T] 
and [2] Panel (a) demonstrates the results for the case when 
the probe electrode is coupled to the interlayer at a; = d/2, 
while the curves, represented in panel (b) are calculated at 
a; = 0. 



Similar to the dc currents, |j°'=''=Ksp) | exhibit quite 
strong non-linearities if Vp is between —(A -I- eV) and 
A + 2eV. The ac spin current also tends to a constant 
value beyond this interval. This limiting value of the cur- 
rent is also a non-monotonous function of the exchange 
field and declines upon increasing h, as for the case of the 
dc spin current. However, unlike the dc spin current, it 
does not contain any contributions, which rise upon in- 
creasing h. The ac spin current is entirely originated from 
the non-equilibrium proximity effect and, consequently, 
is suppressed for the exchange fields, which are consid- 
erably larger than A. For h = 4:.5Tc the ac spin current 
is suppressed by a factor ~ 4 with respect to the case of 
h = 0.9Tc, illustrated in Fig.[3l 

As concerns the limiting behavior of the ac electric cur- 
rent, it is very different from the dc one: it also saturates 
at a constant value instead of to be a linear function of 
Vp. This limiting value of the ac electric current is com- 
parable to or in most part of cases even less than the 
limiting value of the spin current. In particular, in the 
middle of the interlayer (at x = d/2) the limiting value 
of the electric current is negligible (see Fig.[31[a)). There- 
fore, the ac current is highly spin-polarized practically 
for the whole parameter range we consider. We define 
the polarization of the ac current by the following way 



p(ac) 



(17) 



where (...) means averaging over time and P(t) is the 
instantaneous spin polarization of the ac current, defined 
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similar to the polarization of dc current Eq. (|16p: 

\f'''-P{t)\/s'' 



Pit) 



(18) 



The degree of ac current spin polarization, calcu- 
lated according to Eq. (fT7|) . is represented in inserts to 
Figs. [3Ka) andO^b) at a; = d/2 and a; = 0, respectively. 
It is worth to note here that although the ac spin current 
declines upon increasing h, the exchange field also sup- 
presses the ac electric current. For this reason the spin 
polarization of the ac current remains to be high even 
for the case of exchange fields h > A. Another factor, 
which suppresses the proximity effect in the interlayer is 
considerable increase of the voltage V with respect to A. 

As it was already mentioned above, all the most essen- 
tial features of the current flowing through the probe elec- 
trode in the system under consideration (behavior of dc 
spin and ac electric and spin currents, high spin polariza- 
tion of the ac current) originate from the manifestation of 
spin-split proximity minigaps in the non-equilibrium dis- 
tribution function for electrons in the interlayer. In order 
to get insight into qualitative behavior of this distribu- 
tion function one can make use of a balance equation. 
It cannot be applicable for quantitative consideration of 
the discussed problem even in case if the interlayer is 
short enough to justify disregarding the inelastic relax- 
ation processes. Nevertheless, it can help us to clarify the 
qualitative mechanism of the proximity-generated non- 
equilibrium effects, discussed above. Equating flows of 
incoming and outgoing electrons, one obtains the follow- 
ing expression for the spin-up and spin-down distribution 
functions in the interlayer 



N^{e)+NR{e) 



(19) 



Here a =T,i is the electron spin. = /f(£) and 

= fpis — eV) are distribution functions in the left and 
right superconducting leads, fp stands for Fermi distri- 
bution function. N^'^ are local densities of states near 
the interface in left and right superconductors. N'"'^ are 
spin-dependent and exhibit peak features (related to the 
corresponding minigaps in the ferromagnet) for subgap 
energy regions due to proximity effect. Assuming density 
of states in the normal probe electrode to be independent 
on energy and the distribution function to be approxi- 
mately step-like, the current carried by the electrons with 
spin a through the probe electrode is proportional to the 
following expression 



deN^{e,x)f^{e) 



deNAe,x)il~ U{e)) . (20) 



The distribution function and the density of states only 
change considerably within a voltage interval [ei , £2] ■ A 
typical example of the distribution function qualitative 
behavior in the interlayer is represented in Fig. [3I^a) . For 



eV < El all the quasiparticle states are occupied, that is 
fa — I- Otherwise, for eV > £2 all the states are empty 
and fcr = 0. The behavior of the product Na{e)fa{£), 
which is plotted in Fig.dKb), is analogous to that one of 
the distribution function. Therefore, if the potential of 
the probe electrode eVp is taken to be less than £1, then 
the second term in Eq. (|^D|) can be neglected and the 
current carried by spin a electrons is proportional to the 
expression 



(21) 



Npaiei ~ eVp) + / deNaie,x)faie) 



where Npa- is the density of states in the ferromagnet, 
corresponding to large enough quasiparticle energies. It 
is seen from Eq. (pij) that the dc electric current is lin- 
ear function of the applied voltage, as it should be. At 
the same time if one neglects the difference between Np^ 
and Npi , that is only considers the zero order of the pa- 
rameter h/ep, then the dc spin current is only originated 
from the second term, which is independent of Vp. This 
is proximity driven contribution to the current, because 
if the proximity features are absent, this term is spin and 
time independent, as it can be seen from Eq. (|19p . 
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FIG. 4: (a) The qualitative behavior of the distribution func- 
tion faie) in the interlayer in dependence on e. Solid and 
dashed lines are related to spin-up and spin-down electrons, 
respectively, (b) A typical behavior of the product Ncr{e)fcr{e) 
as a function of the quasiparticle energy. 

On the contrary, in the presence of the proximity in- 
duced features the product Na-{£)fcr{e) is different for 
spin-up and spin-down electrons in the interval [£1,62], 
as it is demonstrated in Fig. [DJb) , what leads to the dis- 
cussed above constant limiting behavior of the dc spin 
current. The limiting constant behavior of the ac current 
also originates from this term. The point is that the par- 
ticular shape (especially width) of the minigap features 
is fairly sensitive to the phase difference at SFS junction. 
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As in the nonequlibrium conditions the phase difference 
is driven by the factor 2eVt, all the features generated 
by the minigap in Fig. |4] oscillate in time with the corre- 
sponding frequency, giving rise to ac contribution to the 
current. 

Another interesting property of the ac current flowing 
through the probe electrode is the possibility to obtain 
coherent highly spin-polarized ac currents with a certain 
phase shift between them. The point is that the values 
of the Green's functions cjp'^'^ in the interlayer, which 
determine the current flowing through the probe elec- 
trode, depend on the coordinate x across SFS junction. 
At the same time the frequency of their oscillating part is 
position-independent and controlled by the voltage V ap- 
plied to SFS junction. Therefore, by coupling two probe 
electrode to different locations in the interlayer one can 
implement a source of coherent ac currents having a cer- 
tain phase difference between them. Fig. [5] demonstrates 
the distribution of the ac electric current phases flowing 
through the probe electrode in our system in dependence 
of the coordinate x across the interlayer (vertical axis) 
and the potential of the probe electrode Vp (horizontal 
axis). The phase of the ac current is calculated with 
respect to the phase of the ac electric current, flowing 
through the SFS junction itself, because this reference 
value is independent on Vp and location in the interlayer. 
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FIG. 5: The distribution of the ac electric current phases flow- 
ing through the probe electrode in dependence of the coordi- 
nate X across the interlayer (vertical axis) and the potential 
of the probe electrode Vp (horizontal axis) . 



It is seen in Fig. [5] that the phase shift close to tt can 
be obtained for a wide Vp range by coupling the probe 
electrodes symmetrically with respect to the middle of 
the junction. In principle, phase shifts close to 7r/2 are 
also reachable, although for a narrower ranges of Vp, as 
it is seen in Fig. [51 

The most part of the discussed above results is pro- 
nounced in the regime, when the value of the exchange 
field h (measured in the energy units) is of order of super- 
conducting order parameter A. In ferromagnetic alloys 
like CuNi, which have been intensively used by now for 
experimental investigation of equilibrium properties of 
SFS heterostructures, the exchange field is several times 
larger than A. However, as far as we know, the work 
on the creation of appropriate alloys is in progress now, 
so we believe that this limit can be experimentally real- 
ized in the nearest future. As the F layer is supposed to 
be an alloy, a role of magnetic scattering may be quite 
important^^'S. Therefore, the infiuence of the magnetic 
scattering on the results, obtained in this paper should 
be also investigated. 

In conclusion, in the paper we have theoretically ob- 
tained that voltage-biased SFS junction in the regime of 
the essential proximity effect can be used as a source of 
highly spin-polarized ac current by tunnel coupling the 
interlayer region to the additional normal electrode. This 
current is driven by the non-equilibrium proximity effect 
in the system and vanishes if the proximity features are 
suppressed somehow. The frequency of the ac current is 
controlled by the voltage V applied to SFS junction. In 
addition this system can implement a source of coher- 
ent ac currents having a certain phase difference between 
them by coupling two probe electrode to different loca- 
tions in the interlayer. 
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